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Using mean field theory, we have studied Bose-Fermi mixtures in a one-dimensional optical lat- 
tice in the case of an attractive boson-fermion interaction. We consider that the fermions are in the 
degenerate regime and that the laser intensities are such that quantum coherence across the con- 
densate is ensured. We discuss the effect of the optical lattice on the critical rotational frequency 
for vortex line creation in the Bose-Einstein condensate, as well as how it affects the stability of the 
boson-fermion mixture. A reduction of the critical frequency for nucleating a vortex is observed as 
the strength of the applied laser is increased. The onset of instability of the mixture occurs for a 
sizeably lower number of fermions in the presence of a deep optical lattice. 



I. INTRODUCTION 

One of the most essential differences predicted be- 
tween a superfluid and a normal fluid is the eventual 
appearance of quantized vortices in the former when 
set under rotation Not in vain, since the demon- 
stration of the occurrence of Bose-Einstein condensation 
0- 0- C3| ■ the nucleation of vortices in dilute ultra-cold 
atomic gases in rotating traps has been instrumental to 
elucidate the superfluid properties of these systems and 
has attracted much experimental and theoretical research 
(see e.g. Refs. H S S B H 13 IH IH El and references 
quoted therein) . The prominence of the study of vortices 
in atomic gases transcends the frontiers of the field it- 
self, because topological defects are cornerstones in cross- 
disciplinary areas of physics. In turn, very recently, the 
experimentation with Bose-Einstein condensates in opti- 
cal lattices is stimulating new avenues in the investigation 
of quantum coherence and interference phenomena. 

In the present work we study some properties of a 
slowly rotating Bose-Einstein condensate (BEC) in a 
one-dimensional (ID) optical lattice We are mo- 

tivated by the recent advances achieved in the experi- 
mental manipulation of rotating condensates |15l llq , on 
the one hand, and of condensates loaded into optical lat- 
tices |l7j, on the other hand. Rotating ultracold Bose 
gases in optical lattices are suitable systems to study 
quantum phenomena. In particular, the excitations, sta- 
bility and dynamics of a vortex line in the superfluid 
regime have been recently addressed ^3- Concurrently, 
the physics of degenerate Bose-Fermi mixtures in optical 
lattices is attracting conspicuous attention. Theoretical 
progress is being made in the basic understanding of the 
system's quantum phase diagram and the superfluid to 
Mott-insulator transition [19J , and experimentalists have 
already been able to prepare a degenerate mixture of Rb 
and K atoms in a tight optical lattice pfT ]. 

The aim of this work is twofold. Firstly, we wish to 
address the formation of vortices in a coherent array of 
Bose-Einstein condensates in a ID optical lattice. In do- 
ing this, we will determine the dependence of the thermo- 
dynamic critical angular velocity Q c for vortex formation 
[a, H on the laser intensity. Secondly, we want to inves- 



tigate the effect of trapped fermions on £l c . Attractive 
Bose-Fermi mixtures may experience collapse and thus, 
we will also discuss the impact of the optical lattice on 
the mechanical stability of the mixture. 

The manuscript has been organized as follows. In Sec- 
tion[n]a short description of the mean field model for the 
present problem is given. Section II I II is devoted to the 
discussion of the numerical results for a rotating boson 
condensate confined by magnetic and optical traps. In 
Section ITVl we address the properties of a 87 Rb- 40 K mix- 
ture loaded in an optical lattice. Our concluding remarks 
are laid in Section [V] Finally, a derivation of the virial 
theorem for the mixture in the ID optical lattice is given 
in the Appendix. 

II. THEORY 

We consider a zero temperature mixture made of a 
87 Rb Bose condensate (B) and a degenerate 40 K Fermi 
gas (F). They are confined by the axially-symmetric ex- 
ternal potentials of a harmonic magnetic trap and of a 
stationary periodic optical lattice modulated along the 
z-axis. The lattice is produced by a far detuned laser, 
which hinders the possibility of spontaneuous scattering 
and yields practically equal Rb and K optical potentials. 
The resulting potential for each kind of atom q = B,F 
is 

V q = \ m q (u> 2 q± r 2 + u 2 qz z 2 ) + -j- cos(2nz/d) , (1) 

where m q is the atomic mass and r = \/ x 1 + y 2 is the 
radial variable of cylindrical coordinates. 

The radial and axial frequencies of the harmonic trap 
are taken from a recent experiment fl7| : u>b± = 2tti / b± 
with vb\_ = 74 Hz, and lubz — ^^Bz with Vb z = H Hz 
for 87 Rb, while those for 40 K are a factor (ms/mi?) 1 / 2 ~ 
1.47 larger. The optical potential is determined by its 
period d = 2.7/^m jl7j and depth Vo = sEr, where Er = 
h 2 w 2 /2md 2 — h x 80 Hz is the recoil energy and s is 
a dimcnsionless parameter that provides the intensity of 
the laser beam. 

For 3 x 10 5 atoms of 87 Rb confined in the harmonic 
trap, the condensate is cigar-shaped with a Thomas- 
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Fermi length Ltf = 84 /xm, and a radius Rtf = 6 /xm. 
The superimposed optical potential splits the BEC over 
Ltp/cI ~ 30 wells. We shall be concerned with a range 
of laser intensities that ensures full coherence of the con- 
densate across the whole system The criterion for 
the Mott transition [2l| for the considered parameters 
leads to an estimate of a maximum lattice depth for the 
superfluid regime of V ~ 150 E R ~ h x 12 kHz [l7|. 
Provided that the strength Vq does not surpass largely 
this value, quantum tunneling between consecutive wells 
is still sufficient to retain coherence and one can apply 
the Gross-Pitaevskii theory for the order parameter to 
study the properties of the condensate. 

Within the mean-field approach, neglecting p-wave in- 
teractions, the energy density functional that describes 
the boson-fermion mixture at zero temperature with a 
quantized vortex in the condensate along the z-axis has 
the form El El 
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where k is the circulation quantum number, n B = |^| 2 
is the condensate density, rip is the fermion density, and 
the fermion kinetic energy density has been written in 
the Thomas-Fermi- Weizsacker approximation. The wave 
function of the boson condensate with a vortex line is 
\& = \J n B (r, z) e lK ^, where <f> refers to the azimuthal an- 
gle in the xy-plane. In this work we consider the cases 
k = (non-rotating condensate) and k — 1 (condensate 
with a singly quantized vortex). As is known, vortices 
with more than one quantum of circulation will not exist 
in equilibrium Q. The boson-boson and boson-fermion 
coupling constants g BB and g B F that enter Eq. J2J are 
written in terms of the s-wave scattering lengths a B and 
a B p as g BB — 4ir a B h 2 /m B and g B p = 4:ir a B F?i 2 /m B F, 
respectively, with m B F = 2m B mF / {m B + mp)- 

The variation of £ with respect to n B and n F , un- 
der the constraint of given number of bosons N B and 
fermions Np, yields two coupled Euler-Lagrange equa- 
tions. Namely, a Gross-Pitaevskii equation for bosons 
with a term describing the boson-fermion interaction: 
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and a Thomas-Fermi- Weizsacker equation for fermions 
[22II23I I which is increasingly valid the larger Np is: 
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FIG. 1: (Color online) Contour plots of the 87 Rb condensate 
density in the xz-pl&ne. (a) Condensate in the pure harmonic 
trap, (b) Same as (a) with a quantized vortex along the z- 
axis. (c) Same as (a) under an optical lattice with Vo/h = 
2.5 kHz. (d) Same as (c) with a quantized vortex along the 
z- axis. Distanc es are in units of the oscillator length as± = 
\/h/ (msWBx) = 1.25/xm. In each panel the color scale ranges 
from deep blue at lowest densities to red at highest densities 
[these are 229 (a), 212 (b), 394 (c), and 376 a~ 3 ± (d)]. 



where [i B and fip are the chemical potentials of the bo- 
son and fermion atomic species, respectively. The ground 
state of the boson-fermion mixture is found with k = 0, 
while the vortical state is computed with k = 1. Start- 
ing from randomly sampled densities, we solve Eqs. JSJ 
and J3J by means of the imaginary time method |24l. 
We use the adapted version of the virial theorem 23] to 
check the numerical convergence of the solution (see the 
Appendix) . 



III. BOSON CONDENSATE LOADED INTO AN 
OPTICAL LATTICE 



We first consider a pure boson condensate, with Np = 
3 x 10 5 atoms of 87 Rb as in the experiment of Ref. |17j . 
For the boson-boson scattering leng th we use the value 
a B = 98.98a (la = 0.529 A) 25]. In Fig.HJwe show 
contour plots of the condensate density in the xz-plane, 
n B (x, 0, z) = Y&{x : 0, z)| 2 , for different situations. In (a) 
the condensate is confined only by the harmonic trap. In 
(b) the cigar-shaped condensate of (a) hosts a quantized 
vortex line along the z-axis. In (c) a ID optical lattice 
with Vo/h = 2.5 kHz splits the cigar-shaped condensate 
(a) into an array of multiple disk-like coherent conden- 
sates. 

As shown in Ref. 26], certain aspects of the macro- 
scopic properties and low-energy dynamics of a magnet- 
ically trapped BEC in a tight optical lattice can be un- 
derstood in terms of a renormalized interaction coupling 
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FIG. 2: (Color online) Top: Density profile of the condensate 
along the z-axis ns(0, 0, z) — 1^(0, 0, z)\ 2 , in units of a~^ 3 ± , 
for the indicated values of the lattice depth Vb. Bottom: The 
z-dependence of the combined potential of the harmonic trap 
and the optical lattice, Eq. in units of fkJB±- The hori- 
zontal lines depict the location the chemical potential of the 
condensate boson cloud trapped in that potential. 



constant g* BB > Qbb and of an effective mass m* B > tub 
along the direction of the periodic optical potential. The 
increase of the radial size of the sample that one ob- 
serves in panel (c) with respect to panel (a) is due to the 
increased repulsive effect of the boson-boson interaction 
when the system feels the optical lattice. It originates 
from the fact that the optical confinement produces a lo- 
cal compression of the gas inside each well |2(| . The 
axial size of the condensate also increases appreciably due 
to the combined effect of the repulsive interactions and 
the redistribution of atoms inside the potential wells. 

In panel (d) of Fig. ^ the optical lattice is superim- 
posed to the condensate with the vortex state (b). In a 
bulk superfluid, the size of the vortex core is of the or- 
der of the healing length £ = [^irnBaB]" 1 ^ 2 , where tib 
is the bulk density. This expression holds for inhomoge- 
neous superfluids taking for ns the local boson density 
in the absence of vortices 6]. Thus, the vortex core size 
is larger for the outer sites of the split condensate, where 
the density is smaller. 




FIG. 3: (Color online) Column density ub(x), in units of 
10 3 x a,g 2 ± , for different laser intensities. The column density 
of the condensate when it hosts a vortex line is also plotted. 



In the upper part of Fig. [21 we show the evolution of 
the density profile of the condensate along the z-axis, 
ns(0, 0, z) = 1^(0, 0, z)| 2 , with the laser intensity. In the 
lower part of Fig. we depict the total trapping poten- 
tial felt by the atomic condensate. In the presence of a 
shallow ID optical lattice, the condensate profile starts 
developing small oscillations around the profile at Vb = 
which follow the periodic optical potential, but its spatial 
extension is nearly not affected. When the laser intensity 
increases the system nearly splits into separate conden- 
sates. The effect of the local compression inside each 
potential well of the lattice, as well as the increase of the 
axial size of the system, can be clearly seen in Fig. 

In Fig. 13 we plot the column density i%b(x) — 
J dzriB(x,0, z) for the condensate in the harmonic trap 
alone and when the optical lattice is added. A deep op- 
tical lattice results in a visible reduction of the column 
density at the interior of the condensate and in an ex- 
tended surface. Figure|31also displays ub{x) for the con- 
densate hosting a vortex line. In this case the quantized 
circulation around the vortex line pushes the atoms away 
from it. 

The effect of the modification of the value of the lat- 
tice constant d is illustrated in Fig. 0] In this figure we 
plot the computed density profile ns(0, 0, z) of the boson 
atoms in the condensate for the same parameters as in 
Fig-Elbut for a larger d (d — 6 x 2.7/mi). 

The thermodynamic critical angular velocity fl c for nu- 
cleating a singly quantized vortex is obtained by sub- 
tracting from the vortex state energy E\ in the rotating 
frame the ground-state energy Eq, i.e., Q c = (Ei/Ng — 
Eq/Nb)/^^ |(|, and it provides a lower bound to the 
critical angular velocity 6]. In Fig. we plot f2 c as a 
function of the laser intensity for Nb = 3 x 10 5 and for 
Nb = 5 x 10 4 . For a fixed number of condensate atoms, 
fl c decreases when Vb increases, in agreement with the 
fact that the radial size of the system becomes larger 
in the presence of the optical lattice and the associated 
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FIG. 4: (Color online) Density profile of the condensate along 
the z-axis ns(0,0, z) = ^(O, 0, z)\ 2 , in units of a^, for dif- 
ferent laser intensities when the period of the optical lattice 
potential is increased to d = 6 x 2.7/^m. 
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FIG. 5: (Color online) Critical angular velocity in units of 
ujb± for nucleation of a singly quantized vortex line as a func- 
tion of the laser intensity. 



IV. ATTRACTIVE BOSON-FERMION 
MIXTURE LOADED INTO AN OPTICAL 
LATTICE 



reduction of atoms along the symmetry axis (see the de- 
crease of the column density ns at the origin in Fig. [3] 
with increasing Vo). In turn, the increase of the inertia 
of the condensate along the direction of the laser beam 
due to the larger effective mass m B [2(| also contributes 
to the reduction of the critical angular velocity. 

Assuming that the relative effect remains of the same 
order in experiment, Fig.[3]predicts a sizeable 24% reduc- 
tion of Vt c already with a laser strength Vq/Ji = 5 kHz 
for both of the Nb values. The effect is enhanced the 
shallower the magnetic trap is. A less elongated mag- 
netic trap would favor the appearence of a vortex at a 
slower rotation. For example, the values of f2 c (Vb = 0) 
shown in Fig. [5] would be decreased by a 0.6 factor in a 
spherical trap with vb z — v b±. = 74 Hz. The relative 
reduction of fi c caused by Vo/h = 5 kHz with respect to 
Vb = would still be of 20-23%. As the thermodynamic 
tt c underestimates the actual critical frequency for vor- 
tex creation, dynamical calculations would be needed for 
a precise quantitative estimate of these changes. 

Figure [S] suggests an experiment to test the depen- 
dence of the critical frequency for vortex nucleation on 
the intensity of the optical lattice, since it shows that, 
for a rotating condensate with a given Nb, a vortex line 
should nucleate at angular frequencies much lower than 
^c(Vb = 0) if a co-rotating deep ID optical lattice is su- 
perimposed, whereas otherwise it would not. The exper- 
imental observation of this effect would provide a clear 
signature of the correctness of the physics predicted by 
the Gross-Pitaevskii theory for rotating condensates in 
optical lattices. 



With the purpose of assessing the effects of a trapped 
fermion cloud on the rotating condensate in an opti- 
cal lattice, we consider next a 87 Rb- 40 K mixture. It is 
characterized by a large attractive Bose-Fermi interac- 
tion that assists the sympathetic cooling of the fermionic 
species down to the degenerate regime. There are vari- 
ous values of the interspecies scattering length cibf pub- 
lished, depending on the type of experiment and tech- 
niques used j23, I28I I29I y(| • A recent measurement 
oi cibf at JILA [23 has established a value cibf = 
—250 ± 30 ao, which we adopt here because of its reduced 
experimental uncertainties. Vortex states in 87 Rb- 40 K 
mixtures confined by a harmonic potential have been re- 
cently addressed (23 (although cibf = — 395 ao H^j was 
employed). When a ID optical lattice is switched on, the 
atomic species experience the trapping potential (JIJ. 

We have carried out calculations for a mixture with 
Nb = 3 x 10 5 and Nf = 1.5 x 10 5 atoms, and the same 
parameters as in Fig.Q] We assume that the Fermi com- 
ponent is in the normal — nonsuperfluid — but quantum 
degenerate phase, and consider that it is in a stationary 
state. This situation could be achieved experimentally by 
waiting long enough after the generation of the vortex in 
the condensate to let the drag force between bosons and 
fermions to dissipate. 

We present in Fig. H3 contour plots of the boson 
ns(x,0, z) and fermion np{x, 0, z) densities when the 
condensate hosts a vortex line, for two laser strengths: 
Vo/h = 0.5 (a,b) and 2.5 kHz (c,d). The condensate and 
the fermionic cloud are modulated by the regular pat- 
tern of the optical lattice. The large mutual attraction 
between fermions and bosons makes the effective interac- 
tion between bosons less repulsive than in the pure BEC. 
Furthermore, the boson atoms induce an effective attrac- 
tion between fermions counteracting the Fermi pressure. 
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FIG. 6: (Color online) Contour plots of the condensate (B) 
and fermionic (F) densities in the xz-plane for the 87 Rb- 
40 K mixture, in the combined harmonic and optical lattice 
trap. The condensate hosts a quantized vortex. In (a) and 
(b) Vo/h = 0.5 kHz, whereas in (c) and (d) Vo/h = 2.5 kHz. 
Distances are in units of cib±- In each panel the color scale 
ranges from deep blue at lowest densities to red at highest 
densities [these are 312 (a), 25 (b), 457 (c), and 37 a,g 3 ± (d)]. 



As a result, the density of both species increases in the 
overlapping region, as if they were more strongly confined 
by the external potentials, and the condensate becomes 
more compact in space (compare Figs. [(J: and^l). Still, 
the effect of the Pauli exclusion principle is notorious in 
Fig. [fj]where the 40 K cloud is seen to extend to larger dis- 
tances than the 87 Rb atoms, both axially and radially. 

The presence of a vortex along the z-axis in the con- 
densate component is apparent in panels (a) and (c) of 
Fig. as the boson density vanishes on the vortex line. It 
is interesting to note that the large attractive cibf leads 
to a visible depletion of the fermionic density on the z- 
axis as well, reminiscent of the bosonic vortex core. The 
effect should be directly observable experimentally. It is 
more evident in (d) where bosons and fermions are more 
tightly confined by the optical lattice. 

In the present mixture we find that the presence of 
the quantum degenerate fermionic cloud does not change 
much the value of Q c compared to the pure BEC, nor 
its pattern against the optical lattice strength. Indeed, 
for all the Vo values considered here the critical fre- 
quency for vortex appearence is raised by ~ 10%, due 
to the enhancement of the condensate density at the 
core caused by the fermion atoms. The impact of the 
fermion cloud on f2 c is magnified in a deeper magnetic 
trap. For instance, if the same mixture was set in a trap 
with vbz — 95 Hz and vb±_ — 640 Hz, that preserves 
the aspect ratio vbz/vb± = 0.15 used in Fig. the crit- 
ical frequency f2 c (Vo = 0) would increase by 40% with 
respect to the Np — case. In this squeezing magnetic 
trap, however, the effect of Vq on Q c is negligible (< 2%). 



A peculiar feature of the 87 Rb- 40 K system, that stems 
from the interplay of the moderately repulsive as and the 
strongly attractive cibf, is the existence of a mechanical 
stability limit beyond which the system cannot sustain 
more fermions. At high densities the strong interspecies 
attraction may overcome the Pauli pressure and drive the 
ultracold gas mixture into collapse. The phenomenon has 
been observed experimentally |28j . This motivates us to 
address next the stability of the 87 Rb- 40 K system subject 
to an optical lattice. In fact, collapse is currently the fo- 
cus of experimental attention because of the implications 
the instability has for constraining olbf psl l29l l30l| . and 
also because there are prospects that a radially-squeezed 
87 Rb- 40 K mixture at a density close to collapse must be 
able to form stable bright soliton trains [3l| . 

If we keep Nb fixed at 3 x 10 5 atoms as in the pre- 
vious discussions, with agf = — 250 ao we find that the 
boson-fermion mixture confined in the magnetic trap of 
Fig. would be stable up to virtually arbitrary fermion 
numbers, consistently with the findings of Ref. j^. In- 
creasing u>b±_ prompts the occurrence of collapse. In a 
trap with vg z = 95 Hz and vbi_ = 640 Hz, the 3 x 10 s 
rubidium atoms are able to retain up to a maximum of 
jy-max _ 2 g x 1Q 5 potassium atoms, when the system col- 
lapses. We recall that the effect on Q, c of Np = 1.5 x 10 5 
and Vq for this squeezing trap has been discussed before. 

Now, with vbz = 95 Hz and ub± — 640 Hz, a decrease 
of Np ax can be seen when a ID optical lattice is added. 
It is related to the increase of the effective local compres- 
sion of the atoms caused by the optical trap. Indeed, 
superimposing a ID optical lattice with d — 2.7^m and 
Vo/h = 5 kHz to this strongly confining harmonic trap, 
the 87 Rb- 40 K mixture with Nb = 3x 10 5 can sustain up 
to Aj? ax = 7.5 x 10 4 fermions only, i.e., less than 30% 
of the Vq — value. The effect appears to be accessible 
to verification under present experimental conditions. It 
opens the possibility to study the collapse of a mixture for 
trapped fermion and boson numbers considerably smaller 
than in a pure harmonic trap. 



V. CONCLUSION 

Summarizing, we have determined the thermodynamic 
critical frequency for nucleating a quantized vortex in an 
array of coherent condensates within the Gross-Pitaevskii 
theory. Next we have studied a Bose-Fermi mixture 
trapped in a ID optical lattice to examine the effect on 
the value of Q c of a fermionic cloud in the quantum de- 
generate phase. 87 Rb- 40 K mixtures may collapse and, in 
this regard, it has been shown that a ID optical lattice 
can be an efficient degree of freedom to tune the onset 
of instability. Our analysis of the relative variation of fl c 
calculated by thermodynamic arguments may provide a 
useful information for future experiments. 
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APPENDIX 

In this Appendix we make use of the principle of scale 
invariance to obtain the virial theorem for the energy 
functional of the boson- fermion mixture in an optical lat- 
tice. The virial theorem results from homogeneity prop- 
erties of the kinetic and potential components of the en- 
ergy of the many-body system with respect to a scaling 
transformation that preserves the normalization [3^ . 

A normalized scaled version of the particle density, for 
either bosons or fermions, is 



n\(r) = A 3 ra(Ar) . 



(5) 



where A is an arbitrary scaling parameter. On account 
of ©, a straightforward manipulation shows that the 
various contributions to the energy from the functional 
scale following the rules 



A T 



(6) 



for the kinetic energy (including the boson vortex term 
if present), 



1 

Uhx = -rp Uh 
for the harmonic potential terms, and 

U gX = A 3 U g and U 9bfX = A 3 U g 



(7) 



(8) 



for the intra- and inter-species interaction terms. The 
contribution of the optical lattice transforms as 
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+ ^{Ug BB +U gBF ) 

+ U op t B (Xd) + U op t F (Xd) . 



(10) 



And, finally, the stationarity condition of the energy 
against dilation leads to the virial theorem for the boson- 
fermion mixture: 
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Equation l|ll(l is a useful tool to check the accuracy 
of the numerical solution of the coupled Gross-Pitaevskii 
and Thomas-Fermi- Weizsacker equations. As a typical 
example, in our calculations of the mixture with Nb = 
3 x 10 5 and Np — 1.5 x 10 5 atoms in an optical lattice 
of strength Vq/Ii = 2.5 kHz, from our converged solution 
we find that the rhs of Eq. i|ll|) takes values ~ 9 (~ —80) 
without (with) a vortex line, while Uh is ~ 6 x 10 6 and 



U, 



opt 



is 



— 5 x 10 b . 
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